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Abstract — Sensorless control of Permanent-Magnet Syn- 
chronous Motors at low velocity remains a challenging task. 
A now well-established method consists in injecting a high- 
frequency signal and use the rotor saliency, both geometric 
and magnetic-saturation induced. This paper proposes a clear 
and original analysis based on second-order averaging of how 
to recover the position information from signal injection; 
this analysis blends well with a general model of magnetic 
saturation. It also experimentally demonstrates the relevance 
for position estimation of a simple parametric saturation model 
recently introduced by the authors. 

I. Introduction 

Permanent-Magnet Synchronous Motors (PMSM) are 
widely used in industry. In the so-called "sensorless" mode 
of operation, the rotor position and velocity are not measured 
and the control law must make do with only current measure- 
ments. While sensorless control at medium to high velocities 
is well understood, with many reported control schemes 
and industrial products, sensorless control at low velocity 
remains a challenging task. The reason is that observability 
degenerates at zero velocity, causing a serious problem in 
the necessary rotor position estimation. 

A now well-established method to overcome this problem 
is to add some persistent excitation by injecting a high- 
frequency signal [1] and use the rotor saliency, whether geo- 
metric for Interior Permanent-Magnet machines [2], [3], [4], 
or induced by main flux saturation for Surface Permanent- 
Magnet machines [5]. To get a good position estimation 
under high-load condition it is furthermore necessary to take 
into account cross-saturation [6], [7], [8]. 

The contribution of this paper is twofold: on the one hand 
it proposes a clear and original analysis based on second- 
order averaging of how to recover the position information 
from signal injection; this analysis blends well with a general 
model of magnetic saturation including cross-saturation. On 
the other hand the paper experimentally demonstrates the 
relevance for position estimation of the saturation model 
introduced by the authors [9]. 

The paper is organized as follows: section III] extends 
the saturation model introduced in [9] and studies its first- 



order observability; in section III position estimation by 
signal injection is studied thanks to second-order averaging; 
finally section llV] experimentally demonstrates the relevance 
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of the approach and the necessity of considering saturation 
to correctly estimate the position. 

II. An energy-based model for the saturated 
PMSM 

A. Notations 

In the sequel we denote by Xij := (xi,Xj) T the vector 
made from the real numbers Xi and Xj, where ij can be dq, 
a (3 or 7<5. We also define the matrices 



M„ 



COS /l 

sin/i 



— sm/i 
cos/i 



and /C 



-1 

1 



and we have the useful relation 

^ = K.M, = M.K.. 
dfj, 

B. Energy-based model 

The model of a two-axis PMSM expressed in the syn- 
chronous d — q frame reads 



d<j>dq 

dt 
J duj 

n 2 dt 
dB_ 

dt 



Udq 
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Ri dq - UjK.((j)dq + 4>m) 



idqfc{<i>dq + <j>n 



T_L_ 
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■ U, 



(1) 

(2) 
(3) 



with (j)d q flux linkage due to the current; <f> m :— (A,0) T 
constant flux linkage due to the permanent magnet; Ud q 
impressed voltage and id q stator current; cj and 9 rotor 
(electrical) speed and position; R stator resistance; n number 
of pole pairs; J inertia moment and tt, load torque. The 
physically impressed voltages are u a p := MgUdq while 
the physically measurable currents are i a p := Moidq- The 
current can be expressed in function of the flux linkage 
thanks to a suitable energy function 'ri{4>d,4>q) by 



Idq 



l -dq{<Pdq) 



diU(<j>d, <fiq) 
d2rl((/>d,<Pq) 



(4) 



where dkH denotes the partial derivative w.r.t. the k th 
variable [10], [11]; without loss of generality 'H(0,0) = 0. 
Such a relation between flux linkage and current naturally 
encompasses cross-saturation effects. 

For an unsaturated PMSM this energy function reads 



Hi 



i) 



1 
2Ld 



+ 



2L„ 



where L4 and L q are the motor self-inductances, and we 
recover the usual linear relations 



id = diHi(4> d , 4> q ) 

i q = d 2 Hl(4>d, 4>q) 



L A 



L„ 



Notice the expression for H should respect the symmetry 
of the PMSM w.r.t the direct axis, i.e. 

H(4> d ,-<l> q )=H(fa,<t> q ), (5) 

which is obviously the case for Hi- Indeed (jTJ — <[3]> is left 
unchanged by the transformation 

(u d , U q , 4>d, <f>q, id, lq,U, 9, T L ) -> 

{u d , -Kg, 4> d , -<t>q, id, ~iq, -w, -0, -Tt). 

C. Parametric description of magnetic saturation 

Magnetic saturation can be accounted for by considering 
a more complicated magnetic energy function H, having 
Hi for quadratic part but including also higher-order terms. 
From experiments saturation effects are well captured by 
considering only third- and fourth-order terms, hence 

'rl(cj)d,4> q ) = Hl(<f)d, 4>q) 

3 4 

i=0 i=0 

This is a perturbative model where the higher-order terms 
appear as corrections of the dominant term Hi- The nine 
coefficients cny together with L d , L q are motor dependent. 
But |5]l implies 02,1 = ao,3 = 03,1 = 01,3 = 0, so that the 
energy function eventually reads 

'ri(<pd, 4> q ) = l~Li((j)d, 4> q ) + a 3 . Q (j) d + ai^4>d4> 2 q 

+ a 4 ,0<f>d+ a 2,2 ( t ) d ( l )2 q+ a 0A < t ) t- (® 

From Q and (|6]l the currents are then explicitly given by 

id = -j— + 3a 3fl (f> 2 d + "1,20, + 4a 4 ,o^d + 2a 2 , 2 </>d</>q (7) 

Ld 

iq = -y- + 2a x ^d<i>q + 2a.2,2<t>d<t>q + 4ao.40g, (8) 

which are the so-called flux-current magnetization curves. 

To conclude, the model of the saturated PMSM is given 
by ([l]l-(|3]l and |7]l-(|8]l, with (f) d ,<f) q) u),9 as state variables. 
The magnetic saturation effects are represented by the five 
parameters 0:3,0, 0:1,2, 04,01 02,21 «o,4- 

D. Model with id,i q as state variables 

The model of the PMSM is usually expressed with cur- 
rents as state variables. This can be achieved here by time 
differentiating i dq =Idq(<Pdq), 

didq _ , ^dq 

with -^ 3 - given by ([TJ. Fluxes are then expressed as cf> dq — 
Z dq (idq) by inverting the nonlinear relations |7])-([8]i; rather 



than performing the exact inversion, we can take advantage 
of the fact the coefficients a^j are experimentally small. At 
first order w.r.t. the ay we have <\> d = L d i d + 0(\aij\) and 
4> q = L q i q + 0(\ctij\); plugging these expressions into (j7]i- 
|8]l and neglecting 0(|ai,j| ) terms, we easily find 

4>d = L d (id - 3ot3flL d i d - a lj2 L q i 2 

(9) 



~ ^wL 3 d i d - 


2a 2 , 


2L d L 2 q i d i 2 q 
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Notice the matrix 



Gdd{idq) Gdq{idq)\ nT ( T -l(- \\ 

G dq (i dq ) G q q{i d q) - "^V*, Wh 



with coefficients easily found to be 

1 



Gdd(idq) 



Ld 



6a 3:0 L d id + 12a 4 ,oL d i d + 2a 2y2 L q i 



(10) 



(11) 



-2 -2 



Gdqiidq) = 2a Xi2 L q i q + 4a 2t2 L d idL q i q 

G qq (id q ) = — + 2ai, 2 L d id + 2a 2 . 2 L 2 d i 2 d + l2a ^L 2 i 2 , 



is by construction symmetric; indeed 

nT 11 s _ (dllU{<j> d ,(t>q) d 21 U((j)d,<j>q) 
U± d qK<Pdq) ~ y dl2U{4>d: g ) d22 H(cbd, 4>q) 

and <9i 2 "H = d 2 \H. Therefore the inductance matrix 



Ldd(idq) L dq {i dq )\ (G d d{i d q) Gdq(idq) 

L dq {idq) L qq (i dq )J \Gdq{idq) G qq {idq) 

is also symmetric, though this is not always acknowledged 
in saturation models encountered in the literature. 

E. First-order observability 

To emphasize position estimation is not easy at low speed, 
we study the observability of ([TJ — <[3j, augmented by 



dr L 
dt 



= 0, 



(12) 



around a permanent trajectory defined by 



= u dq - Ri dq - uJC(4> dq + 4> m ) 

= ^dq^&dq + <t>m) ~ ~ 
d9 

— = u. 
dt 

Notice such a permanent trajectory is not a steady state point 
unless UJ := since 6 hence u a p and i a p are time-varying 
(<f> dq ,i dq ,u dq ,uJ,TL are on the other hand constant). Also 
the rationale for augmenting the model with ( fT2| ) is that it is 
usually also desired to estimate an unknown load torque. 
From i a p = Mgi dq we then get 

5i a p = SM e idq + Mg6i dq = M- s (K,i dq 89 + 6i dq ), (13) 



and similarly for u a p. The linearization of <JXJ — <(3j and ( p"2| ) 
around a permanent trajectory is then 



dS(f>dq 



dt 

J dSui 3 
rf~dt 



Mj(Su a p - R6i a/ 3) - K.((/) dq + 4> m )8ui 
+ Lu(((j) dq + <f> m )66 - IC64> dq ) 



dS6 

~dt 

d&TL 

dt 



= ^KpM^(^ dq + <t> m ) 

~ ^dqi&dq + <t>m)S6 - K.S(j) dq ) — 

= 8u> 



= 0. 



where we have used u dq — Ri dq = U)K{<j) dq + 4> m )- 
On the other hand time differentiating 8i a g yields 

d8i aj3 dMg / dS9 dSi dq 



LLlICSiap + MgUCidqSw + Dl dq ( 



>dq) 



, d8(t> dq 
dt 



where we have used 5i dq — DI dq ((j> d )Scj) dq . Therefore 



lDi dq g, dq )]-Hii d -^ 

= [DT dq (ci> dq )]~ 1 lCi dq 8uj 



dq 



+ LC(Si a p,Su afj ) 



dqh -" vuv "" ' dt 
= JC<P{4> dq )Suj + <S>((j) dq )Uj59 + LC(8i a0 ,Su a0 ) 

where 



*( 



'dq 



,J dq 



K,[Dl dq {cf> dq )]- x Kl t 



'-dq(Pdq) 



and LC(5i a fj,Su a fj) is some matrix linear combination of 
5i a s and 8u a p. Similarly 



J dSoj 3 - St l 



■LC(6i a p,6u a p). (14) 



Now &(<f> dq ) is non-zero in any reasonable situation, hence 
&(4> dq ) and K<&((/) dq ) are independent vectors. If uj ^ it 
is thus clear that 89 and Sio can be expressed in function 
°f ia/3, ~ d j^~ an d u a p\ as a consequence 8tl is a function 



d<)i 



d-Si 



of i a p, "^, a d l t 2 P and u Q/3 , — ^ by ( fl4) ; finally 



?<((, 



is by (jT3j a function of 8i a p. In other words the linearized 
system is observable by 8i a p. 

If u = only <5w can be recovered from i Q/ g, ^ n(i 
and u a p\ as a consequence only <&{4> dq )88 + ^- can be 
recovered from (14) , and nothing new is gained by further 
differentiation. In other words the linearized system is not 
observable, as pointed out in [10]. 

Estimating the rotor position at low speed is thus inher- 
ently difficult. Yet it is doable when the motor exhibits some 
saliency, geometric or saturation-induced, with the help of 
permanent excitation such as signal injection. Indeed a more 
thorough analysis along the lines of [12], [13] would reveal 
the system is in that case observable in the nonlinear sense 
provide the motor exhibits some saliency. 



III. Position estimation by high frequency 

VOLTAGE INJECTION 
A. Signal injection and averaging 

A general sensorless control law can be expressed as 



U al 3 : 


= Me c u~fS 


(15) 


d9 c 
dt 


= u c 


(16) 


dr\ 
~dl. 


= a(M g J 7 s,9 c ,r],t) 


(17) 


UJ C ■ 


= n c (M g J lS ,9 c ,r],t) 


(18) 


u 7 s : 


= Ujs(M0 c irys,O c ,r],t), 


(19) 



where the measured currents i a p = MgJ-s are use d to 
compute M 7 a, uj c and the evolution of the internal (vector) 
variable r\ of the controller; 9 C and lo c are known by design. 
It will be convenient to write the system equations ([T])-([3]) 
in the 7 — 8 frame defined by x 1 s '■— Me-e c x dq , which gives 

d(j}js 

dt 
J dui 

n?~dt 
d9 

dt 
where from (f7|-([8| currents and fluxes are related by 



U 7 s - Ri 1 S - Ucfc&fS - uJCMg-e^m, 


(20) 


3 Tj 

-^/C(<?f> 7 ,5 + M g - 9c 4> m ) 


(21) 


u, 


(22) 



i 7 «5 = M e ^gld q {M g _ 0a (j) lS ). 



(23) 



To estimate the position we will superimpose on some 
desirable control law (jT9j a fast-varying pulsating voltage, 



U lS ^U 7 s(MgJ 7 s,0 c ,7],t) +Uy S f(Q,t), 



(24) 



where / is a 27r-periodic function with zero mean and 
u^s could like U^s depend on Mg c i 1 $,9 ci ri,t (though it is 
always taken constant in the sequel). The constant pulsation 
Q, is chosen "large", so that j(Vtt) can be seen as a "fast" 
oscillation; typically Q, := 2ir x 500 rad/s in the experiments 
in section HVl 

If we apply this modified control law to (p0|-(p2|i, then it 
can be shown the solution of the closed loop system is 



UjS 



b lS = ^ s + ^F(nt) + 0( — ) (2b) 



1 



* = * + o(- 2 ) 
e c = e c + o(~) 



(26) 
(27) 
(28) 
(29) 



where F is the primitive of / with zero mean (F clearly 
has the same period as /); (cj) lS ,uj,9,9 c ,fj) is the "slowly- 



varying" component of (^ s , u, 9, 9 C , rj), i.e. satisfies g ^ £(j) ._ I f T R ^ ^ g)ds = Q 

fas - - - - ° 

= « 7 <j — Ri~ ( s — u c JC<f) yS — LolCMg_g^ <fi m We have set 



dt 
J duj 

r^~dt 
d9 
dt 
dh 
dt 
drj 

~df 






7<5 



M- a 



n 



= a(M- { 



ig Us,9 c ,r],t), 



where 



(30) 



i jS = M- B _- e T dq (M- 6 _QJ> lS ) 

u^s = U~fS (MgJ 7 s,9 c , f], t) . 

Notice this slowly-varying system is exactly the same 
as (p0]i-(p2]i acted upon by the unmodified control law (fl"5|l- 
([19j. In other words adding signal injection: 

• has a very small effect of order O ( ^ ) on the mechan- 
ical variables 9, to and the controller variables 9 c ,rj 

• has a small effect of order O(^) on the flux <j) 1 $; this 
effect will be used in the next section to extract the 
position information from the measured currents. 

The proof relies on a direct application of second-order 
averaging of differential equations, see [14] section 2.9.1 
and for the slow-time dependance section 3.3. Indeed setting 
£ '■= jl and x := (</> 7 ,5, w, 9, 9 C , rj), p0|)-(|22|) acted upon by 
the modified control law (fl5]>— (p~8]> and ( 24 1 is in the so-called 
standard form for averaging (with slow-time dependance) 

— = eh{x,ea, a) := e(f 1 (x,ea) + f 1 (x,ea)f(a)), 

with fi T-periodic w.r.t. its third variable and e as a small 
parameter. Therefore its solution can be approximated as 

x(a) = z(a) + e(ui(z(a, ea, a)) + 0(e 2 ), 

where z(a) is the solution of 

dz 

da 
and 

I " T 
T 



5l(y,£cr) 
vi(y,scr,a) 



egx(z,ea) +e 2 g 2 {z,ea) 



f 1 (y,ea,s)ds = f x {y,e(j) 



(h(y,£v,s) - gi{y,ea))ds 



h(y,sa) / f(s)ds 



1 f T 
ui(y,ea,a) :=vi(y,ea,a)- — / v 1 {y,ea,s)ds 

= f 1 (y,ea)F(a) 
K 2 (y,ea,a) 
■= difi{y,ea,a)m{y,ea,a) - dim{y, ea, a)gi{y,ea) 

dF 2 (a) 



[fi,fi](y,£°)F(a) + -difi{y,sa)f 1 (y,ea) 



We have set 

[7iJi](y, z°) ~difi(y, £<j)h(y, ea)-d 1 f 1 (y, eajj^y, ea) 

and F(a) := J^ f{s)ds — ^ J Q J^ f(s)dsda, i.e. F is the 
(of course T-periodic) primitive of / with zero mean. 

Translating back to the original variables eventually yields 
the desired result (f25l— d29b. 



B. Position estimation 

We now express the effect of signal injection on the 
currents: plugging (|2"5j) into (f23|) we have 



1-yS 



Ah a 






+0{ i)^7* 



|^) + <^ 



1 



= i T j+i Ti F(ftt) + oy 



(31) 

where we have used ( [30] l and performed a first-order expan- 
sion to get 



V ■■= M-e-e c m d q [,Mj_- e ^ s )Mi - n 



rp U-yS 

fT 



Mg_g^DXdq 



{^( M Le c 



o<^) M leM- (32) 



We will see in the next section how to recover i 7 s and i 7 s 
from the measured currents i 7 g. Therefore ( |32] ) gives two 
(redundant) relations relating the unknown angle 9 to the 
known variables 9 c ,i c i q ,i 1 s,u c [ q , provided the matrix 

rT 



S(y,,i 7 s) 



:= M^Dl, 



dq i^dq ( M J V' 



M, 



effectively depends on its first argument y. This "saliency 
condition" is what is needed to ensure nonlinear observabil- 
ity. The explicit expression for S(y,,i 7 s) is obtained thanks 
to ( fTT) . In the case of an unsaturated magnetic circuit this 
matrix boils down to 



S{n,i lS ) = M^ 



h 



M, 



L d -L 



2LdL q 



1 L d +L 
Ld — L q 
Ld+L q 



- cos 2/1 
i 

sin 2/i 



1 



Ld — L q 
Ld+L q 

Ld — 



Ld+L, 



sin2jU 
cos 2/z 



da 



and does not depend on i-ys', notice this matrix does not 
depend on /i for an unsaturated machine with no geometric 
saliency. Notice also ([32} defines in that case two solutions 
on ] — 7r, 7r] for the angle 9 since S(y, i^s) is actually a 
function of 2/i; in the saturated case there is generically only 
one solution, except for some particular values of i 1 s. 

There are several ways to extract the rotor angle informa- 
tion from ( (32) , especially for real-time use inside a feedback 
law. In this paper we just want to demonstrate the validity 
of pZ) and we will be content with directly solving it through 
a nonlinear least square problem; in other words we estimate 
the rotor position as 

~ 2 

= c -\- arg min z 7(5 - S(fi, ijs)^r 



(33) 



TABLE I 
Rated and saturation parameters of the test motor. 



Rated power 


1500 W 


L d 


7.9 mH 


Rated current I n (peak) 


5.19 A 


L q 


8.2 mH 


Rated speed 


3000 rpm 


&3,oL^I„ 


0.0551 


Rated torque 


6.06 Nm 


ai,2L d L q I n 


0.0545 


n 


5 


"4,0irf^ 


0.0170 


R 


2.1 O 


d2,2L d L 2 q ll 


0.0249 


A 


0.155 W b 


a 0,4 L Pn 


0.0067 



C. Current demodulation 

To estimate the position information using e.g. ( (33] l it is 
necessary to extract the low- and high-frequency components 
i 7 s and i-ys from the measured current i 1 g. Since by pTj ) 
i 7 s{t) ~ i-ysif) +i~/s(t)F(nt) with i 7 g and i 1 s by construc- 
tion nearly constant on one period of F, we may write 

rt 



i~/s(t) 



i-yS(t) 



1 



t-T 



i- f s(s)ds 



f t _ T ijs(s)F(fls)ds 



-10C- 



-120 
-14C- 



£ 



where T := 4?, Indeed as i* 1 is 2tt -periodic with zero mean, 

/ i 7 s(s)ds ~ i 7 g(t) / ds + i^s(t) / F(fls)ds 
Jt-T Jt-T Jt-T 

= Ti lS (t) 
/ i lS {s)F{VLs)ds w* 7 *(t) / F(ns)ds 

Jt-T Jt-T 

+ V(*) / F 2 (ns)ds 

Jt-T 

r T 

= i-y S (t) / F 2 (fls)ds. 
Jo 

IV. Experimental Results 

A. Experimental setup 

The methodology developed in the paper was tested on 
a surface-mounted PMSM with parameters listed in table [I] 
The magnetic parameters were estimated using the procedure 
of [9, section III]. Notice this motor has little geometric 
saliency (Ld ~ L q ) hence the saturation-induced saliency 
is paramount to estimate the rotor position. 

The experimental setup consists of an industrial inverter 
(400 V DC bus, 4 kHz PWM frequency), an incremental 
encoder, a dSpace fast prototyping system with 3 boards 
(DS1005, DS5202 and EV1048). The measurements are 
sampled also at 4 kHz, and synchronized with the PWM 
frequency. The load torque is created by a 4 kW DC motor. 

B. Validation of the rotor position estimation procedure 

1) Control law: Since the goal is only to test the validity 
of the angle estimation procedure, a very simple V/f open- 
loop (i.e. fi c and U 7 s do not depend on i 1 $) control law 
is used for (fi~5]l-(fT9|i; a fast-varying (ft :— 2tt x 500 rad/s) 
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Fig. 1. Long test under various conditions: (a) measured 8 — 8 C , estimated 
8 — 9 C with and without saturation model; (b) measured speed ui, reference 

,rd 



speed ui c ; (c) load torque t^; (d) voltages u rc i 



square voltage with constant amplitude is added in accor- 
dance with d24l), resulting in 



LO c (t) 



d6c 
dt 
u lS = «$(<) + u c (t)4> m + U 7 <5/(Ot). 

Here U) c {t) is the motor speed reference; u r A(t) is a filtered 
piece-wise constant vector compensating the resistive voltage 
drop in order to maintain the torque level and the motor 
stability; u^s '■— (u, 0) T with u :— 15 V. 

The rotor position is then estimated according to ( [33) 1. 

2) Long test under various conditions (Fig. IT): Speed 
and torque are slowly changed over a period of about two 
minutes; the speed remains between ±2% of the rated speed 



and the torque varies from 0% to 150% of the rated toque. 
This represents typical operation conditions at low speed. 

When the saturation model is used the agreement between 
the estimated position 9 and the measured position 6 is very 
good, with an error always smaller than a few (electrical) 
degrees. By contrast the estimated error without using the 
saturation model (i.e. with all the magnetic saturation pa- 
rameters otij taken to zero) can reach up to 70° (electrical). 
This demonstrates the importance of considering an adequate 
saturation model. 

3) Very slow speed reversal (Fig. |2J: The speed is very 
slowly (in about 20 seconds) changed from —0.2% to +0.2% 
of the rated speed at 150% of the rated torque. This is 
a very demanding test since the motor always remains in 
the first-order unobservability region, moreover under high 
load. Once again the estimated angle closely agrees with the 
measured angle. 

V. Conclusion 

We have presented a new procedure based on signal injec- 
tion for estimating the rotor angle of a PMSM at low speed, 
with an original analysis based on second-order averaging. 
This is not an easy problem in view of the observability 
degenaracy at zero speed. The method is general in the sense 
it can accommodate virtually any control law, saturation 
model, and form of injected signal. The relevance of the 
method and the importance of using an adequate magnetic 
saturation model has been experimentally demonstrated on a 
surface-mounted PMSM with little geometric saliency. 
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